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Abstract 

Since 2008 the VLBI group at the Institute of Geodesy and Geophysics at TU Vienna has focused 
on the development of a new VLBI data analysis software called VieVS (Vienna VLBI Software). One 
part of the program, currently under development, is a unit for parameter estimation in so-called global 
solutions, where the connection of the single sessions is done by stacking at the normal equation level. 

We can determine time independent geodynamical parameters such as Love and Shida numbers of the 
solid Earth tides. Apart from the estimation of the constant nominal values of Love and Shida numbers 
for the second degree of the tidal potential, it is possible to determine frequency dependent values in 
the diurnal band together with the resonance frequency of Free Core Nutation. In this paper we show 
first results obtained from the 24-hour IVS R1 and R4 sessions. 


1. Introduction 

The Vienna VLBI Software VieVS (Boehm et al., 2009 [1]) is a data analysis software, which 
is developed by the VLBI group at the Vienna University of Technology. The software is written 
in Matlab with the plan to make it compatible with equivalent non-commercial software, e.g., 
Octave. One of our main goals is to develop an easy-to-handle software, which will be attractive 
to our students and thus allow us to easily involve more students in VLBI research. In the 
modeling part of the software (vie_mod) the theoretical VLBI observable (group delay) is computed 
following the most recent IERS Conventions (McCarthy and Petit, 2004 [7]) and IVS standards 
like the treatment of thermal deformation of VLBI radio telescopes (Nothnagel, 2009 [9]). The 
time varying parameters (such as Earth orientation or troposphere parameters) are estimated using 
piecewise linear offsets at integer fractions of integer hours with a least-squares algorithm to easily 
allow comparison and combination with estimates from other space geodetic techniques. The 
VieVS software is extended by an extra program unit for parameter estimation in so-called global 
solutions, where the connection of the single sessions is done by stacking the normal equations. 
In this paper we focus on the determination of geodynamical parameters of the solid Earth tides: 
frequency dependent Love and Shida numbers in the diurnal tidal band together with the Free 
Core Nutation period. 

2. Global Solution in VieVS 

Because of the limited computer memory capacity it is essential to keep the equation system 
small. In the VLBI technique there are parameters in the observation equations which cannot be 
fixed to a priori values, even if we are not interested in them (e.g., clock parameters). Therefore, 
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a reduction algorithm is used which is based on a division of the normal equation system into 
two parts. The first part contains parameters which we want to estimate and the second part the 
parameters which will be reduced: 

( N u N 12 W dx! \ ( h \ 

V N21 N 22 )'[dx 2 ) [b 2 ) ’ W 

where N = A T PA and b = A 1 PI. The reduction of dx 2 is done by executing the matrix 
operation in Equation (1) and introducing the second equation into the first one: 

(-N11 — N^N^Nk) ■ dx 1 = b\ — N\ 2 N 22 b 2 N re duc ' dx\ = b re d uc ■ (2) 

Stacking is used for combining normal equation systems if a parameter is contained in at least 
two normal equation systems and only one value in the resulting combined system should be 
estimated. For a combined solution of the identical parameters (dx), the normal matrices (N re d uc ) 
and the right hand side vectors ( b re duc ) from n single sessions have to be summed up: 

NreDUC — Al re. due. A T NreducJi T ••• T -Areduc-n > (3) 

bREDUC b rf:r iuc 1 T b re d uc 2 T ••• T b re duc-n • (4) 

The final solution is obtained by an inversion of the normal matrix: 

dx i = Areduc • bREDUC ■ (5) 

3. Love Numbers and Free Core Nutation Period 

Solid Earth deformation arises from the variations in the Earth’s gravitational field caused by 
the Moon and the Sun. Love and Shida numbers ( h , l ) are dimensionless parameters, which reflect 
the amount by which the surface of the Earth responds to the tide-generating potential (V tul ). 
Considering the Earth being spherical, non-rotating, elastic, and isotropic — i.e., the most basic 
model — the tidal displacement (u) at a given latitude and longitude (ip, A) in the topocentric 
system is described by equations (6)— (8), where g stands for gravity acceleration and n is the 
degree of the tide-generating potential. 
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When a more precise model of the Earth is considered — i.e., ellipticity, rotation, and the elastic 
mantle/fluid core boundary is taken into account — the relationship between the tide- generating 
potential and the displacement becomes more complicated. The displacement vector is composed 
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of deformational parts ( 5u ) of specific harmonic degree and order, and specific frequency inside 
the band (Wahr, 1981 [14]). 

5u R(f) = sin ( 2 7>) sin(0/ + A) , (9) 

where Hf is the Cartwright- Tayler amplitude of the tidal term, <5/i2i(/) stands for the difference 
of the frequency dependent Love number from the constant second degree Love number I 12 , and 9 f 
is a tidal harmonic argument. The strong frequency dependence of Love numbers in the diurnal 
band arises from the resonance behavior of the Earth, which is caused by the presence of the 
fluid core. The rotational axis is slightly inclined w.r.t. the axis of rotation of the mantle. In this 
situation, forces arise at the elliptical core/mantle boundary, which try to realign the two axes. 
In the terrestrial reference frame, this phenomenon is seen as a diurnal motion, in the celestial 
frame it appears as a retrograde motion of the celestial pole with a period of approximately 430 
days and is designated as Free Core Nutation (FCN). Because this effect is a free motion with 
time-varying excitation and damping, resulting in a variable amplitude and phase, an FCN model 
is not included in the recent precession-nutation model adopted by the International Astronomical 
Union (IAU 2006/2000A). It means that after taking into account the precession- nutation model 
a quasi-periodic unmodeled motion of the celestial pole in the celestial frame at the 0. 1-0.3 mas 
level still remains in the measured data (e.g., Lambert, 2007 [6]). 

3.1. Estimation of Love Numbers and FCN in VieVS 

Basically, three approaches can be used to determine the period of the Free Core Nutation 
motion. The first option is an indirect estimation from the analysis of the celestial pole offsets 
obtained from VLBI measurements (e.g., Herring et ah, 1986 [4], Vondrak and Ron, 2006 [13]). 
The FCN effects are also seen in gravity data, where the first determination of the resonant period 
was done by Neuberg et al. (1987) [8] and was followed by, for instance, Sato et al. (1994, 2004) [11], 
[12] or Ducarme et al. (2007) [2]. In the VieVS software we implemented the ability to determine 
the FCN period directly from the analysis of the observed solid Earth tidal displacements of the 
VLBI antennas, as it was introduced by Haas and Schuh (1996) [3]. For the estimation of the FCN 
period a resonance formula for the Love numbers in the diurnal band is used, which was published 
by Wahr (1981) [14]: 


h-2i(iOT) = h2i{iooi) + has • — — • ( 10 ) 

^fcn — wr 

^ 21 (^ 01 ) is the Love number of the 01 tidal wave, whose frequency looi is used as reference 
which is sufficiently distinct from the resonance frequency of the FCN ujfqn- This formula is 
applied in the least-squares adjustment as a condition equation constraining the estimates of 
diurnal Love numbers at the tidal waves T, by means of fitting them to the resonance curve, for 
which also the resonance strength factor has was estimated. Because the equation is not linear 
w.r.t. the ujpcNi iterations have to be carried out. 

4. Data and Results 

For this paper we used all IVS R1 and R4 sessions, i.e., 24-hour observing sessions which 
have been performed two times per week since January 2002 using global VLBI networks. Before 
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the analysis we removed outliers from the observations and excluded sessions with an a posteriori 
variance of unit weight higher than 1.5. Clock parameters, zenith wet delays, troposphere gradients, 
and Earth orientation parameters were reduced from the normal equation system and estimated 
implicitly. The station coordinates and velocities of the 18 VLBI radio telescopes together with 
the radio source coordinates were fixed to their a priori catalog values, ITRF2005 and ICRF2, 
respectively. We solved for Love numbers of ten diurnal tidal waves using partial derivatives from 
Equation (9) and the resonance parameters (i.e. , FCN period with the resonance strength factor 
lifts) were estimated after five iterations from the condition equation (10). 

Our estimates of the Love numbers (Figure 1, 
Table 1) are in very good agreement with the 
recently adopted IERS values. Larger differences 
can be seen mainly for the tidal waves in the 
vicinity of the resonance period (phil, psil) and 
for the tides with lower amplitudes (e.g., NOl, pil, 
Jl). For the resonance period we obtained 413 ± 
20 sidereal days in the celestial frame. The large 
standard deviation of the FCN period might be 
due to using only eight years of data and the high 
correlation with the resonance strength factor, as 
was pointed out in Haas and Schuh (1996) [3]. 
An apparent time-variation of the FCN period 
was discussed in Roosbeek et al. (1999) [10] or 
Hinderer et al. (2000) [5]. 


Frequency dependent Love numbers in diurnal band 



Figure 1. Estimates of the diurnal Love num- 
bers /121 and the resonance curve. 


Table 1. Estimates of frequency dependent Love numbers for ten diurnal tidal waves, FCN period, and 
resonance strength factor compared with the IERS values [7] and with the results from Haas and Schuh 
(1996) [3]. 


Frequency [°/h] 

Tide 

IERS values [7] 

Love number 
this paper 

h 21 

Haas and Schuh (1996) [3] 

13.3987 

Q1 

0.6033 

0.6033 ± 0.0065 

0.560 ± 0.012 

13.9430 

01 

0.6026 

0.6026 ± 0.0017 

0.606 ± 0.002 

14.4967 

NOl 

0.6004 

0.6005 ± 0.0158 

0.435 ± 0.030 

14.9179 

pil 

0.5882 

0.5895 ± 0.0456 

0.623 ± 0.090 

14.9589 

PI 

0.5823 

0.5842 ± 0.0027 

0.574 ± 0.005 

15.0411 

K1 

0.5261 

0.5385 ± 0.0049 

0.496 ± 0.002 

15.0821 

psil 

1.0439 

1.2023 ± 0.1587 

-0.136 ± 0.228 

15.1232 

phil 

0.6623 

0.6590 ± 0.0626 

0.702 ± 0.121 

15.5126 

TT1 

0.6113 

0.6104 ± 0.0830 

0.934 ± 0.158 

15.5854 

Jl 

0.6105 

0.6096 ± 0.0159 

0.538 ± 0.031 


FCN period 

432 sid. days 

413 ± 20 sid. days 

426.3 ± 20.3 sid. days 

llRS 

- 

-0.00217 ± 0.00015 

-0.00162 ± 0.00016 
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5. Summary 

The VLBI group at Vienna University of Technology is developing a new data analysis software, 
called VieVS (Vienna VLBI Software). The first version of the software has been released and can 
be downloaded by registered users. In this paper we focused on a new module for global adjustment 
of more sessions, which is implemented in the software. Besides estimation of a new terrestrial 
reference frame, celestial reference frame, and Earth orientation parameters, the program unit 
allows the determination of geophysical parameters such as frequency dependent Love and Shida 
numbers and the Free Core Nutation period from solid Earth tidal deformations. 
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